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SUMMARY 

A procedure for interpolating between specified points of a curve or sur¬ 
face is described. The method guarantees slope continuity at all junctions. 
A surface panel divided into p x q contiguous patches is completely specihed 
by the coordinates of (p -|- 1) x (g -|- 1) points. Each individual patch, how¬ 
ever, depends parametrically on the coordinates of 16 points, allowing shape 
flexibility and global conformity. 
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I. Interpolation Between points on a Space Cnrve 


We assume that a space curve is approximately defined by a sequence 
of points {A, 5, C,each of which are vectors, i.e., A = Ax,Ay,Az in a 
Cartesian reference frame. We propose an interpolation scheme that defines 
the curve between each pair of adjacent points so that the curve is as smooth 
as possible. Consider four adjacent points: D,E,F,G. 


Figure 1 


The distance t along the chord between E and F defines a scalar variable 
that we use to parametrize a curve c(t) between E and F. c(t) will be a blend 
of two parabolas. 

The three points V, E, F, define a parabola p{r) as follows. Let r be the 
distance along the chord between D and F. Let u be the distance perpen¬ 
dicular to r in the plane defined by D,E,F. The parabola, 

u = ar{d — r), (1) 

has an axis perpendicular to the line along r. d = distance between D and 
F. We choose a so that this parabola passes through E. Then p(r) can be 
taken to be the parabola, Eq.(l). 
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Similarly, the three points E, F, G, dehne a parabola q{s) as follows. Let 
s be the distance along the chord between E and G. Let v be the distance 
perpendicular to s in the plane dehned by E, F, G. The parabola, 

v = j3s{e-s), (2) 

has an axis perpendicular to the line along s. e = distance between E and 
G. We choose f3 so that this parabola passes through F. Then q{s) can be 
taken to be the parabola, Eq.(2). If E,F, and G happen to be collinear, 
Eq.(2) will be n = 0, a straight line. 

We now define c{t) as a blend of p{r) and q{s) as follows: 

c{t) = [1 - {t/to)]p{r) + (3) 

to = distance between E and F. Consequently, the two blending functions 
[the coefficients of p{r) and ^(s)] vary linearly between 0 and 1. Eq.(3) is not 
completely specihed until we dehne a relation between t and r and between t 
and s. This can be done only by dropping perpendiculars to the lines along 
r and s respectively; so r = r(t) and s = s(t) are dehned by this geometric 
operation. 

In like manner we construct a curve Ci(tj) between each adjacent pair 
of points. It is easy to prove that where two Ci(ti) connect, their slopes 
are equal; so the entire curve will be continuous and smooth. To see this, 
consider Eq.(3) rewritten as follows. 

c{t) = p{t) + {t/to)[q{t) - p{t)]. (4) 

The slope at point E is 

{dc/dt)E= {dp/dt)E + (t/to)E[dq/dt - dp/dt]E + {l/to){q- p)e- (5) 


3 



The second term on the right hand side is zero because t = 0 at E; and the 
third term on the right hand side is zero because q = p at E. Consequently, 

{dc/dt)E = {dp/dt)E- ( 6 ) 

The slope of c(t) at E equals the slope of the parabola p{r), which passes 
through D, E, and F. This same parabola determines the slope of c at E for 
the curve between D and E by an identical argument. Consequently slope 
continuity is assured. 

Blending of two parabolas is possible only if the interval is an interior 
one. If the curve starts at point A, 


Figure 2 


then interpolation between A and B should be by the single parabola, defined 
as above, through A, B, and C. Accordingly, in order to specify carefully the 
shape of the curve near its endpoint, the distance between A and B should 
be smaller than that for interior intervals. There is no requirement that the 
points be equally, or nearly equally spaced. Obviously point density should 
be higher where the curvature is higher. 

If an interpolated segment of the curve does not behave according to some 
a priori intent, an extra point can be inserted to provide the required con¬ 
straint. Aside from endpoints discussed above, it should not be necessary to 
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specify for special treatment interior points at inflections or cusps, provided 
the latter are not intended to be infinitely sharp. Precision is reproducing 
such features can be assured merely by specifying extra points sufficiently 
near to these characteristic features. 

The curve c(t), Eq.(3) or (4), is of cubic order in the Cartesian coordi¬ 
nates. Consequently it can adequately represent an inflection that should 
occur within an interval. It should be apparent that the cubic function ap¬ 
propriate to the interval EF does not pass through points D and G. (This is 
an important difference in comparison to simple £t by cubic polynomials.) 
The manner of construction (as a blend of two parabolas) guarantees that 
spurious wiggles will not be introduced, as frequently happens when simple 
cubes are forced to pass through four points of a curve. 

II. Interpolation Between Points on a Space Surface 

We assume that a space surface is defined by a net of points Aij, such as 
intersections of lines in the figure below. 


Figure 3 
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Connecting lines between adjacent points of the net can be constructed by the 
blending algorithm described in Sec. I. Clearly a good surface interpolation 
for the patch EFGH must take into account the global shape implied by 
the adjacent points, L, M, N, ...V,W. Consider the patch EFGH. A set of 
coordinates to which points on the patch are related can be dehned, say x 
and y. X = constant (between 0 and 1) dehnes a space curve through the four 
points a,b,c,d which, in the notation of Sec. I, are points t = xto of curves 
MN, EE, GH, VU respectively. The line i, between b and c, is then defined 
by the blending algorithm. Similarly, a line m is dehned hj y = constant 
(between 0 and 1). Surface point z{x,y) can be taken to be 


z{x,y) = i{y), 

(7) 

z{x, y) = rh{x). 

(8) 


These two surfaces will differ slightly since i and m will not in general in¬ 
tersect. They do coincide, however, on all network lines. If computational 
labor is not a factor, one may take the average of (7) and (8): 

1 - 1 

z{x,y) =-i{y) +-m{x). (9) 

An advantage of these surface interpolation schemes is that all surface slopes 
will be continuous at boundaries between patches. The surface shape of each 
patch will depend, as it should, on the behavior of the surface surrounding 

it. 

Parabolic blending along i, m lines is possible only for interior patches. 
Interpolation for patches at an edge will employ a single parabola for i or m, 
like the end interval of Sec. I. Gorner patches will necessitate single parabolas 
for both i and m lines. Accordingly good edge dehnition requires a narrow 
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spacing for the network lines adjacent to edges of the surface, as shown in 
the hgure. 

An entire surface panel consisting of pq patches is specihed by the coor¬ 
dinates of (p -|- l)(g + 1) points, slightly more than one point per patch on 
the average. The shape of a single (interior) patch, however, depends on 16 
points. This signihes both flexibility and global relationship. 

The optimum arrangement of network points will naturally depend on 
the shape of the panel that is to be represented. If the panel has a ridge, for 
example, the network should be arranged so that a network line passes along 
the crest. Closely spaced network lines (approximately parallel to the crest) 
should also be included to provide adequate dehnition. 

III. Analysis 

In order to carry out the procedures described in the preceding sections 
it is necessary to have explicit formulas for the parabolic functions and co¬ 
ordinates. We need to this only for one case. 


Figure 4 


J is the point along DF obtained by dropping a perpendicular from E. If J 
is -|- x{F — D), then: 

{E-[3+ x{F - D)]} ■ {F - D) = 0. (10) 
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It follows that 


x = {E-D)-{F- D)/£ 


( 11 ) 


where d= \F — D\. The equation of a point p on the parabola p(r) is 

p{r) = D + {r/d){F — D) + ar{d — r){E — J). (12) 

The coefficient a is determined by requiring p{xd) = E. The last term of 
(12) must satisfy, 

axd{d — xd){E — J) = E — J. (13) 

So, 

a = l/[d‘^x{l — x)]. (14) 

The equation of the parabola is now completely specihed: Eq.(12) together 
with values for a and x, Eqs.(14) and (11). 

The only problem that remains is to hnd the relation between t and r. 
From the hgure it follows that 

r = xd + tcos9. (15) 

We have, 

cos e = {F-E)-{F- D)/dto (16) 

where to = \F — E\. Consequently, Eq.(15) is specihed completely, so Eq.(12) 
can be written as a function of t. 

In the alternative event that t is along the chord from D to E, so that 
to = \E — D\, we have instead of (15), 


r = t{E-D)- {F-D)/dto. 


( 17 ) 
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